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Abstract 

Simple states, such as isobaric analog states or giant resonances, embedded into con- 
tinuum are typical for mesoscopic many-body quantum systems. Due to the coupling 
to compound states in the same energy range, a simple mode acquires a damping width 
("internal" dynamics). When studied experimentally with the aid of various reactions, 
such states reveal enhanced cross sections in specific channels at corresponding resonance 
energies ("external" dynamics which include direct decay of a simple mode and decays 
of intrinsic compound states through their own channels). We consider the interplay 
between internal and external dynamics using a general formalism of the effective nonher- 
mitian hamiltonian and looking at the situation both from "inside" (strength functions 
and spreading widths) and from "outside" (5-matrix, cross sections and delay times). 
The restoration of isospin purity and disappearance of the collective strength of giant 
resonances at high excitation energy are discussed as important particular manifestations 
of this complex interplay. 
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1 Introduction 



Dynamical features of open mesoscopic quantum systems are characterized by the presence of 
"simple" (single-particle and collective) excitations, "complicated" (chaotic) intrinsic motion 
involving many degrees of freedom, and irreversible decay into continuum. The coexistence 
and interplay of these phenomena is the important aspect of all processes including excitation 
and deexcitation of the system. One of the questions of primary interest in nuclear physics, 
especially for future experiments with radioactive nuclear beams, is that of existence and purity 
of simple modes of nuclear excitation embedded into continuum. Similar problems arise in 
atomic and molecular physics, physics of atomic clusters and mesoscopic solid state devices. 

During the last decade, a number of related phenomena were discovered in this area of 
nuclear physics, see for example ^ |2|]. Saturation of the spreading width of the giant dipole 
resonance (GDR) in hot nuclei 0, |], "disappearance" of the collective strength of the GDR 
at high excitation energy ||, and existence and relatively narrow width of the double GDR 
|7|, ||, || are J us t few bright examples. In the physics of isobaric analog states (IAS), one can 
mention the evidence for existence of the so-called "broad poles" [|10||, very weak fluctuations 



of the spreading widths of the IAS throughout the periodic table |TT], [12|, [L3| , and restoration 
of isospin purity at high excitation energy ||. 

In such problems, one always deals with a simple excitation mixed with the dense back- 
ground of complicated states. The simple excitation is associated with a specific signal. It 
can be a quantum number which singles out the state in the ocean of surrounding states as 
it happens in the IAS case. As a rule, such a state is relatively pure with respect to this 
label when looked at in the entrance channel. The isospin purity is violated by the internal 
mixing |TT| when, due to the high background level density, the statistical enhancement of 
perturbations becomes extremely important, similar to the well known enhancement of weak 
interactions observed in parity nonconservation |14| . The individuality of a " simple" mode can 
also be referred to its specific structure, for example in the case of a giant collective vibration, 
whose coherence makes the state very different from the background. Such a special state is 
characterized by a large multipole moment which provides a strong collective gamma decay 0. 
In all cases, the manifestations of the simple mode in specific reaction channels are intertangled 
with the chaotic mixing inside the system. 

As a result of the mixing, the simple mode is fragmented over exact stationary states which 
form the fine structure of the spectrum. Being averaged over the unresolved fine structure, 
the excitation function is related to the strength distribution of the original "label" smoothly 
depending on excitation energy. More detailed statistical analysis of observed fluctuations, 
assuming generic correlations of energies and strengths for the invisible underlying states, is 
capable |15j of extracting their characteristics. In general, the strength functions and reaction 
cross sections represent two sides of the process, internal and external, and the relation between 
them is far from trivial. Thus, the strength distribution may or may not coincide with the width 
distribution seen in the reactions and decays Jl6| . 



The well known formalism W7\ of the strength function proceeds as if the states under 
consideration were stable. However, all excited states, strictly speaking, have a finite lifetime 
and therefore belong to the continuum spectrum. The level widths of the resonances in the 



continuum [18|, [L9| are governed by the interaction which is in general different from that 
forming the discrete spectrum inside the system. The effects of intrinsic mixing and coupling 
to and via continuum have to be considered simultaneously. 

Below we formulate a consistent quantum-mechanical approach which fully accounts for the 
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interplay of internal damping and decay and contains, as particular cases, the 
of the collective strength of the giant resonance |lj explained by kinetic arguments in 
the restoration of isospin purity at high excitation energy [0 in accordance with the old idea by 



'disappearance 




and 



Morinaga [EJJ and Wilkinson | 22"|| . We discuss the general properties of the strength function of 
a simple mode embedded into continuum in its relation to what is observed in reactions. Our 
consideration, being intentionally schematic and less specific than in the well known review 
paper fll|, is in many aspects complementary because of its generality and the simultaneous 



treatment of internal and external aspects of the problem. 



2 Effective hamiltonian 



We use the effective non-hermitian hamiltonian |19], El] in order to take into account internal and 



external interactions on equal footing. The intrinsic structure at high level density produces the 
set ("background") of the basis intrinsic states \n), n = 1, . . . N, where A is supposed to be 
large. The simple state |0) is located in the the same range of energy. All N + 1 states have the 
same values of exact integrals of motion such as total angular momentum. We assume that the 
basis states are characterized also by quantities like isospin or parity which are approximate 
integrals of motion. The isospin mixing which is one of the subjects of our application is 
introduced explicitly by the off-diagonal elements of the hamiltonian. Parity nonconservation 
due to weak interactions can be another example of an approximate conservation law which 
can be included in a similar manner. 

The effective hamiltonian in (N + l)-dimensional space is the operator 

n = H- % -w (i) 

containing two (real and symmetric for a time reversal invariant system) matrices H and W 
which describe internal and external coupling, respectively. 

The antihermitian part W has a special structure JL9|, |TT| EJ]] being originated by the on-shell 
decays into open channels c = 1,2, ... ,k, 

W = AA T W nn , = £ A c n A c n ,. (2) 

c 

Here we introduced the (N + 1) x k matrices A = {A c n } of real transition amplitudes which are 
proportional to the matrix elements of the full original hermitian hamiltonian which connect 
intrinsic and channel subspaces of total Hilbert space. 

The hermitian part H consists of the unperturbed energy e of the simple state |0), the 
internal N x N hamiltonian h describing the background states \n), and the coupling between 
the simple and complicated states. The real coupling matrix elements H 0n = H n0 = V n ,n > 1, 
form an A- dimensional column vector V. The spectrum h n of the eigenvalues of h is supposed 
to be very dense. Along with the similarity of generic complicated wave functions ||25|| , this 
justifies the statistical approach. 



The effective hamiltonian 7i can be studied with standard matrix methods [Ej , |26| . Its 
diagonalization gives complex eigenvalues 

£j = E J - t -r J , j = 0,...,A, (3) 
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and the quasistationary eigenstates \j) with a pure exponential decay law ~ exp(— iSjt). The 
construction of the effective hamiltonian guarantees the unitarity of the scattering matrix, see 
below Sect. 5.2. 



3 Standard model of the strength function 

The description of the mixing of stable internal states, which forms complicated stationary 
superpositions and spreads the srength of original simple states, is well known |fT7|] . With the 
antihermitian part W omitted, the intrinsic propagation within the closed system is described 
by the Green function G(E) of the hermitian part of the hamiltonian, 

G(E) = ^ . (4) 



The eigenvalues of the intrinsic hamiltonian H are given [17j by the (N + 1) poles E = e a of 
the Green function (f|). They are the roots of the secular equation 

G Q0 \E) F(E) = E - e - £ l^V = «• ( 5 ) 

n>l h - h n 



Each eigenfunction \a) of if carries a fraction 



r = |c, 



o|2 
I 



(6) 



of the collective strength determined by the weight of the corresponding component Cq in the 
expansion over the basis states, 

N 

H=^|0) + £C<». (7) 

n=l 

The smooth strength function of the simple excitation is defined in terms of the average local 
level spacing D of background states, 

P (e) = [r/D{e)\e a =, (8) 

It is normalized according to J2 a f a = I deP (e) = 1. 

The formal solution (|6]) requires the knowledge of statistical properties of the background 
spectrum h n and coupling matrix elements V n . The simplest ansatz used in the standard 
model [17[ assumes a roughly equidistant dense spectrum of h n and interaction intensities 
V% uncorrelated with energies h n and slightly fluctuating around their mean value (V 2 ). For 
convenience of the reader, we collect the results of the uniform model in the Appendix, along 
with a brief discussion. 

At (V^ 2 ) > D 2 , the strength function of this uniform model has the Breit-Wigner shape 

P ° (e) = ^( e -e ) 2 + (ri)V4- (9) 
with the spreading width given by the golden rule, 

(V 2 ) 

= r s = 2tt^Z. (io) 
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The standard model just described is valid [^7], |2D| if T s does not exceed the energy range AE 
of coupling strength V% (AE is defined by the spread of the doorway states which provide the 
gates for the further mixing of the state |0)). This is expected to be a good approximation for 
the IAS with the typical spreading width < 100 keV. In the case of giant resonances T s ~ AE 
and the uniform model should be corrected [30]. However the difference influences mainly the 
shape of the wings of the strength function which is of minor importance for our purpose; here 
we use the uniform model for definiteness. 



4 Simple state embedded into continuum 
4.1 Formulation of the problem 

Now we take into account the openness of the system. The simple state |0) is open to the direct 
decay (channels c which display specific signatures of the simple mode, for example collective 7- 
radiation from the giant reasonance or pure isospin of the IAS). Due to the intrinsic coupling to 
compound states, the simple state also acquires access to many "evaporation channels" labeled 
by the superscript e; partial widths depend on the distribution of strength of the simple mode 
carried by specific compound states. 

When applied to the IAS with isospin T>, we have to consider the surrounding background 
states \n) which belong mainly to the isospin T — T> — 1. The isospin mixing occurs mostly 
through intrinsic interaction [|TT| so that the decay channels for the decoupled simple mode and 
evaporation channels for compound states carry different isospins. In many cases, the effects 
we are interested in can be studied using one direct channel which will be labeled as c = 0. 
Then we have in the hamiltonian (||) the amplitudes A® = v /7o, where 70 is the "natural" width 
of the simple state, and A e n1 n > 1. All A e n are assumed to be of the same order of magnitude. 

At low energies (for example, for neutron resonances), only few decay channels are open and 
the narrow compound states do not overlap. Their widths 7„ = J2e(^n) 2 are small compared to 
their mean energy spacing D. As energy and level density increase, we pass the region of strong 
coupling via continuum where the width collectivization occurs and broad "Dicke resonances" 
[ P3L P^ , form the contribution of direct processes. The situation changes again when many 
uncorrelated decay channels are open, and the off-diagonal elements (0) of the antihermitian 
part W of the effective hamiltonian are averaged out. Then the states |0) and \n) simply acquire 
finite widths. These unstable states are coupled through the hermitian interaction V and this 
is what bridges the gap between the intrinsic strength function and its manifestation in the 
resonance reactions. 

To describe the open compound states, we introduce the N x N Green function 

z — h + (i/2)w 

N complex poles z = e u of g(z) determine energies and evaporation widths of compound reso- 
nances still decoupled from the simple mode. In eq.(|ll]), w stands for the (N x N) submatrix 
of W, eq.(|), which acts in the compound subspace and describes the evaporation together 
with the interaction between the compound states through common decay channels. The latter 
is characterized by the off-diagonal matrix elements of w. Being the sums of uncorrellated 
contributions of many evaporation channels, k ^> 1, these elements, due to mutual cancella- 
tions, are small in comparison with the diagonal elements, \w nn >/w 

nn 

l/Vk (see H). The 
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corresponding corrections are of order of 7 2 1 ,//cD 2 where ^ ev is the typical evaporation width. 
We will neglect them below assuming ^ ev <C y/kD. Under this condition, partial decay widths 
of the compound states to specific evaporation channels are small, lev/k <C D. 

The complex energies of compound resonances in this approximation are equal to e„ = 
h v — (i/2)j ev ,is = 1,2, ...,N, supposing on the statistical grounds that the fluctuations of 
widths of compound states are weak since the number k of evaporation channels is large. The 
simple state has its own complex energy e = e — (2/2)70 where 70 is the direct decay width. 

Let us now switch on an interaction between the simple and compound states through the 
hermitian coupling operator V. The mixing proceeds in competition with the decays of intrinsic 
states, both via direct and evaporation channels. Whence, we need to generalize the standard 
procedure of determination of the strength function, Sect. 3, for the decaying system. In our 
schematic although quite generic model, it could be done exactly. 



4.2 Decay widths in the presence of intrinsic damping 

The diagonalization of the total non-hermitian hamiltonian (|I|) leads to N + 1 complex eigen- 
values (H) which are the roots z = Sj of the secular equation (compare to (||)) 

F(z)=Z-e O -V T g(z)V=0, (12) 

or, in the explicit form, 

^-^-E7Tr = - ( 13 ) 

The interaction amplitudes V u , which couple the unstable simple state \0) with complicated 
(and decaying as well) intrinsic states are still real in the approximation taken above (we 
neglected the off-diagonal part of the continuum coupling w). 

Similar to (0), the quasistationary eigenstates \j) can be represented as superpositions of 
decoupled unstable states \0), . . . , \v) 

\j) = cm + T.c>). (14) 

V 

The fraction fi = |Cq| 2 of the strength of the simple state \0) carried by the quasistationary 
state \j) is equal, as in eq.(0), to 

~P = xhj' LJ = v V(^-M^)v. (is) 

With Sj = Ej — (i/2)Tj, the loops V can be written as 

V = V - V » , 9 = — ImVl/ 2 — L^. (16) 

v 1*^3 L 3 lev v t»j c u 



Using the secular equation (|i3[) we arrive at a very simple expression 



V = lLJ± f (17) 

J- j lev 



leading to the individual strengths (|T|) 



ft = Tj lev . (18) 

70 - lev 
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In other words, the resulting width of the quasistationary state \ j) can be found from simple 
probabilistic arguments, 

T 3 = l0 ~P + lev {l-p). (19) 

The direct decay width is distributed over all quasistationary states according to their fractions 
of the strength of the original simple state. It is easy to check the normalization of the weights 

E/ i = ^ r (Er J -(^ + i)7«)=i (20) 

where the last step follows from the invariance of the imaginary part of the trace of the hamil- 
tonian (|l|), J2j Fj = 7o + Nj ev . We have to notice that the probabilistic interpretation emerges 
here as a result of a strict quantum-mechanical calculation, with no ensemble averaging or 
transition to a kinetic description. 



4.3 "Broad pole" 

Explicit expressions for the properties of the strength function, including the spreading width 
along with the decay widths into continuum can be obtained if the average characteristics of 
the intrinsic spectrum and of the coupling matrix elements are specified. 

In the uniform model |T7| used earlier for the stable states, eq.(|T3|) gives a pair of coupled 



equations for the real and imaginary parts of the complex energy (^), see Appendix, 

E 3 =e, + -T s fif , r 3 - = 7o-i y fj . 21 
2 1 + xjyj 1 + x]y] 



Here T s stands for the standard spreading width (0), and notations 



^ = cot (tt-^J , % . = tanh ^_^_^j (22 ) 

are introduced. 

At moderate values of the interaction V, the simple state keeps an appreciable fraction of the 
collective strength and preserves its individuality, see eq. (E4f) below. Such a state was called 



the broad pole in |10| . The problem of IAS can serve as a typical example. The unperturbed 



analog state \0) arises at the energy eo carrying almost pure isospin T>. Its direct decay width 



7o is much larger than the evaporation width 7 e „ of background states with isospin T< |19 
The isospin-violating interaction V mixes these states. 

Assuming that the resulting width r = Tj =0 satisfies the condition (r — j ev ) >I>we have 



from ( |22|) y « 1, so that eqs. (plf) give for the complex root corresponding to the broad pole 

E « e , T = 7o - r s . (23) 



The second expression reads r = T* — T^- in the notations chosen in [ 10| . The state T> can be 



observed only if it decays before mixing, 70 > T s . The collective strength (|1^) carried by the 
broad pole is then 

f° = ^° ~ ^ ev ~ ^ s — 1 _ ^ s (24) 

which remains of order of unity as long as T noticeably exceeds 7 e „. This formula extends 
to the case of unstable compound states the measure introduced in |L0j, |3l[] of the "purity of 
analog spin" of the broad pole. On the other hand, the typical values fi for j 7^ are small. 
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4.4 General strength function 



The energy dependence of the strengths fll8p is hidden in the secular equations (pi). Exclusion 
of Xj leads after simple algebra to the general equation for the strength function which depends 
only on the absolute value I70 — j ev \, 



P 



D 



2vr|7o - 7 e 



In 



eo) 2 + \ 



r s + |7o-7e*|(l-/ j 



or, for small /. 



P 



D 



{E 3 



In 



^o) 2 + \ 



|7o-7e«|(l-/ J ') 



2 ■ 



{Ej 


-e ) 2 + i(r s + | 7 o 


-7e,|) 2 


(Ej 


-e ) 2 + i(r s -| 7 o 


-7e,|) 2 



(25) 



(26) 



2tt|7o - 7 e 

Substituting summation over j by integration over energy, one can easily check that this 
distribution is normalized as 

r. 



dE H(F 



7o - 7e 



1, 

l-r. 



T s < 70 - 7e 



T s > |To - 7e 



(27) 



■7e«| (compare with eq. (|23|)) of the simple 



In the upper case, the contribution / — l — l s / 17 

state (|23"D is lost in the integral. Indeed, the small factor in front of the logarithm in eq. (|25|) 
is in this case compensated due to the small value of the denominator of the expression under 
the logarithm so that eq. (|26|) is not valid for this special state. On the other hand, when 
the increasing mixing rate characterized by the spreading width prevails upon the influence of 
direct decays, the simple mode fully dissolves in the sea of compound states. 



Except for an exponentially narrow domain of parameters around the point r s 
the width (FWHM) V of the distribution (|26| ) is determined by 



|7d -lev\, 



ir 



(70 - lev f 



{21 



The tails of the strength function, E 3> (T s + 17 — J e v\), are universal and given by the standard 



model, p w (D / '2n)T s / ' E 2 . In the limits T s > 
the Breit-Wigner distribution 



7 eu | or Y s < |7o - 7e„|, eq.p|) reduces to 



D 
2^ 



r. 



< 



2 D 



(^■-e ) 2 + r2/4 - 7T r s 





7o - 


lev 




\l0 - lev 


(Ej - e 


) 2 + 


70 " 


- lev 


2/4 



2 

< - 



(29) 



respectively. Near the point T s 



: |7o 
D 
7rT a 



le 



7T (70 - 7e^) 2 ' 

(p6|) is invalid and eq.(^) gives 



ln- 



27Tl\ 



D 



In In 



27rr fl 
D 



+ ... 



(30) 



The strength f° is still larger than all / J for j ^ but this cannot influence the normalization 
fl2"7p. Fig. 1 illustrates the relation between the exact expression for the strength function, 
eq.(^), the approximation (^), which is invalid in the center of the spectrum, and more crude 
approximations (p9|). 

The strength function gives an average description of the fragmentation of individual simple 
configurations in the intrinsic space. In the next section we study the problem as it is seen in 
continuum properties. 
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5 Scattering characteristics 



5.1 Scattering matrix 

Up to now we concentrated on the "inside" view of a simple unstable mode mixed with com- 
plicated fine structure states. The "outside" world was present as a reservoir for irreversible 
decay through numerous open channels. Now we take a glimpse of the same system from the 
viewpoint of reaction amplitudes and cross sections where only asymptotic states are observed. 
The scattering matrix S = {S cc } at energy E can be written as |19| 



S(E) = s 1 ' 2 {l-if{E)}s 1/ \ (31) 

f(E) = A T Q(E)A. (32) 

Here s includes the potential scattering as well as channel coupling and direct reactions in the 
continuum. Those effects being unrelated to intrinsic dynamics are irrelevant for our purpose, 
and s(E) can be considered as a diagonal matrix with phase shift elements exp(2i5 c ) smoothly 
depending on E. The Green function in eq.(|3~2"D, 



S{z) = j-^, (33) 

describes the propagation governed by the total hamiltonian (|I|). It differs from the intrinsic 
Green function by the antihermitian part of the effective hamiltonian. Both Q{z) and the 
scattering matrix (|3"T| ) have poles at the complex energies (|3|). 

It is a straightforward exercise to establish, with the aid of the factorized structure (2) of 
the antihermitian part W, the relation 

Q(E) = G(E) - -G(E)A A T G(E) (34) 



between the two Green functions ([33|) and (f|). The i?-matrix in eq . (|3^) is familiar from nuclear 
reaction theory P^ ], 

R(E) = A T G(E)A. (35) 

It describes the propagation inside the closed system between two acts of coupling to the 
continuum; the poles of R(E) correspond to the energies e a of intrinsic states with the mixing 
V fully accounted for. The reaction matrix T(E) of eq.(^) is similar to fl3"5]) but includes 
all intermediate couplings to the continuum. Finally, for the scattering matrix ([H], ^2|) the 
substitution fl3~3| ) gives 

f{E)= ^1 , S(E) = s^ 1 - m ^ E) ^. (36) 
l + (i/2)R{EY l + (i/2)R(E) 

5.2 Scattering wave function, delay time and unitarity 

The scattering wave function \^ C E ) with the incident wave in the channel c at energy E can be 
presented by the superposition of intrinsic, |n), and continuum channel, \c; E), components, 

r°° i 

l^> = E^(£)N+E / ,dEY c (E,E')\c';E') (37) 
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where E c ' is the threshold energy in the channel c'. Recall that the decay amplitudes A c n are 
the matrix elements of the total original hamiltonian between the states |n) and \c;E). By a 



direct substitution of ( |37|) into the Schr0dinger equation, we find |19| the N x k matrix h(E) 
of the intrinsic components b c n as 

b(E) = G(E)As 1/2 . (38) 

The diagonal elements of the k x k matrix b'(E)b(E) determine the norm of the internal 
part of the wave function initiated in the channel c at energy E. Therefore this matrix should 
characterize the fraction of delay time in this reaction due to intrinsic resonances. Indeed, the 



Smith's time delay matrix is defined [R3] as 



f(E) = -i&(E)^f. (39) 
Taking into account only the resonance energy dependence via the i?-matrix in eq. (|36"D , we find 

f res (E) = -S-W U s 1 ' 2 . (40) 

l-(i/2)R(E) dE \ + {i/2)R{E) 

In the same resonance approximation one can neglect the energy dependence of amplitudes A 
to get from (j35|) 

(dR/dE) res = —A T G 2 {E)A. (41) 
Using the relation (|3"4| ) between the Green functions G and Q, we obtain 

f res = r 1 / 2 A T g\E)g(E)As 1 / 2 . (42) 



Thus, the time delay matrix (|40| ) coincides with the intrinsic norm matrix found from ([38]), 

b f b = free- (43) 



The total Green function ( |3~3"|) describes the propagation in the open system and, therefore, the 
delay time as well. 

We can now define the normalized probability p c n (E) to find the system in the intrinsic state 
\n) in the "elastic" reaction c — * c, 

Pn{E) = -^-K(E)\ 2 , J2Pn(E) = l. (44) 

res\ ) n 

The probability Pq{E) characterizes the weight of the simple state |0) in the channel c. In the 
problem of the IAS this quantity measures the isospin purity in a given channel. 



The full scattering matrix (|31~| , |3~2"D is unitary provided the potential scattering matrix s 
is unitary. It follows from the fact that the decay amplitudes A in the entrance and exit 
channels of eq.(|32|) are the same which appear in all intermediate processes described by the 
total propagator Q(E) with the aid of the effective hamiltonian (|l]||). 

The unitarity condition SS' = S^S = 1 gives for the reaction matrix ( p2|) 



f tf = i(f - f f ) (45) 
which can be transformed, with the help of (|2|) and (p8|), into 

s l/2 b\E)Wb(E)r 1/2 = i{f{E) - T\E)}. (46) 
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6 A simple case: Stable background states 



The simplest situation corresponds to the stable background states with no direct access to 
open channels, 7et) — > 0, when the intrinsic evolution for the reaction in the channel c starts 
and ends at the simple state. The background states are involved by the internal coupling only 
at the intermediate stages of the reaction. Calculating the diagonal element of the resonance 
time delay matrix (32) we obtain for the probability (H), 



pl{E) = f(E) = [l + E (F V \ J '' = W/dE]-\ (47) 

This is nothing but the continuous generalization of the strengths f a = \Cq\ 2 defined above by 
eq.(^]) in discrete points e a of the intrinsic energy spectrum, f a = f(E = e a ). Since the intrinsic 
states are coupled to continuum through the state |0) and the probabilities p c n are normalized, 
eq.([P]), the decay (or population) partial widths 7 q do not appear in (f4"T|). If several direct 
decay channels c are open, the energy behavior fl4"7|) is identical for all of them being determined 
by intrinsic dynamics only. 

The probability fl4"T|) vanishes at energies E — h n of the unperturbed background states 
which are located intermittently with the actual energies e a . In the vicinity of h n the compli- 
cated states dominate the intrinsic part of the scattering wave function. 



Another, though equivalent to (|47|) , representation of the time delay in terms of the complex 



energies (|3|) of quasistationary states can be derived from ( f42|) 



r Z (E) = -2|MW(S)1 = f E ~ E _ ^ + m (48) 

The delay times for different channels c are proportional to the corresponding partial widths 
of the state \0) and have the identical energy dependence determined by the complex energy 
spectrum of intrinsic unstable states. 

The representation (^) is useful when the vicinity of the broad pole (|23| ) is considered. It 
follows from (^) that the contribution of this pole is a smooth function of energy superimposed 
onto the picket fence of the 5-like peaks with the average value proportional to the Weisskopf 
recurrence time tt/D for a long-lived wave packet. At the energy E = Eq, the time delay in a 
channel c due to excitation of the broad pole is equal to 47Q/70IV On the other hand, one gets 

l^™i 2 = 7o C |^ooTO| 2 ^4 7o 7( 7o ) 2 (49) 

since the energy E is very close to the unperturbed energy of the state |0). Therefore the 
probability maximum is determined by the fraction of the total width 7o of the original mode 
which still resides at the broad pole, 

Po(^o) = — = • (50) 

7o 7o 

in agreement with (^4|) taken at ^ ev = 0. 

One should have in mind that the distribution (|47|) wildly fluctuates on the fine structure 
energy scale. With the energy resolution worse than the level spacing D, one sees only a 
smooth behavior coinciding with that of the strength function Pq(E), eq.(|8|). It is quite natural 
because here the intrinsic mixing is the only source for the spreading of the strength, or for 
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isospin impurity in the case of IAS. An average magnitude of the probability to find the original 
isospin can be easily estimated in the standard model with the uniform background. Eq.(|47|) 
gives here (1^ = T s ) 

. . sin 2 (Ett /D) . . 

sm {E7r/D) + [7Tl l /2D) 

or, after averaging over fine structure, and taking 1^ ^> D, 

1 D 



ME) = -Tr r (52) 

This natural estimate (inverse number of fine structure states within the spreading width) 
coincides with that used by von Brentano |10| . 



7 Mixing with open compound states 
7.1 Purity of a simple state 

The situation changes in the realistic case with many open evaporation channels. Strong 
fluctuations of the probability Po(E) are smeared out since the compound poles are displaced 
to the complex energy plane even with no coupling to the simple mode. This probability 
remains considerable in a finite vicinity of the point eo ensuring a noticeable isospin purity of 
the internal part of the scattering wave function in this region. 

If the simple mode and the compund states have no common decay channels, the nonzero 
decay amplitudes are Aq = y^TJ (consider for simplicity a single direct decay channel) and A e n . 
The reaction amplitudes are equal to 



T CC (E) = lo g 00 (E), T ce (E) = v/to £ g Qv {E)Al (53) 

V 

where now Q^(E) = T{E) (see ([12])) whereas 

Gov(E) = -^Goo(E). (54) 
h, — e v 

The delay time in the elastic process, according to ( fl2j) and (J38[) , is given by 

r c r c es {E) = 7o|goo(£)| 2 fl + E | E ^ ja ) (55) 



= l0 \g 00 (E)\ 2 [l + L(E)] 



where the loop L(E) is the analog of L 3 , eq. flTSp , taken at the running real energy E rather 
than at the complex energy Sj. Therefore we find instead of ([17]) 



1 



l + L(E) 



V: 



p Q (E) = f(E) = - ■ = 1 + • (56) 



The function f(E) extends the strength function (O) of the quasistationary states to a running 
real energy E (compare with the similar correspondence between the functions (@) and ([17]) in 
the case of stable compound states). Note that, by definition (fi"4"P, the resonance envelope \Goo\ 2 
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is divided out of normalized probabilities po{E) which behave uniformly within the spreading 
width. 

The loop function ( |56|) can be calculated similar to fllTf) . Under the same assumptions, it 
is equal to 

7e„ 1 + 2 V 



^H-^T^- (57) 



where, instead of u22|), we now have 



x = cot(vr— ), y = tanh(--^). (5i 



For a small evaporation width, 7 e „ <C D, the expression (|57|) reduces to 



L(E) = ^(l+x 2 ). (59) 

The results in this limit do not depend on the evaporation width at all and therefore coincide 
with those following from (pET|). In particular, the weight of the simple state in the intrinsic 
part of the scattering wave function is in average of order of D/Y s <C 1. 

As level density and number of open channels increase, the ratio , y ev /D rapidly grows to- 
gether with the argument of y, eq. fl5"8|) . One has a fast transition to the limit of the overlapping 
background states when y pa 1 and L(E) — > Y s /^ ev . The probability fl56|) in this case is 
noticeably greater than in (|52"D, 

lev ~ 1 s 1 s 

The fluctuations disappear, and the simple state preserves its individuality in the intrinsic wave 
function across the whole region of the giant or analog resonance. This behavior is demonstrated 
in Fig. 2. 

The purity of the intrinsic part becomes perfect when 7 et) ^> r s ; the depletion of admixed 
states of the opposite isospin occurs faster than their population. This gives a microscopic 
justification of the isospin purity at high excitation energy predicted in |H, 22] and recently 
observed experimentally |2j. At the same conditions, the fraction of the simple mode carried 
by a generic compound state, 

1-/=— Xr"> ( 61 ) 

lev + 1 s 

is small. This means that the compound processes have no time to explore the presence of the 
exceptional simple state. 

The equivalent result was formulated in terms of the kinetic balance between the processes of 
decay and mixing in [ 20] where the mechanism for the disappearance of the collective strength of 
the GDR at high energies was suggested. The authors showed that the probability of excitation 
of a collective mode in an initially heated nucleus is equal, using our notations, to T s / (j ev + r s ) 
and therefore diminishes as r s /7 e „, when the temperature exceeds a critical value determined 
by the condition ^ ev ~ T s . Complementary to the somewhat qualitative kinetic arguments of 
[p0[ , here the analogous conclusion follows from a full quantum-mechanical consideration. 

We need to mention parenthetically that such statements assume the saturation of the 
spreading width pa T s as a function of temperature. The absence of a considerable depen- 
dence on excitation energy is well known for the IAS [11] , |l2l , |13 , [Slf . The saturation of the 



intrinsic spreading width presumably takes place for the GDR as well [0, [3], |], General theo- 
retical arguments in favor of such a saturation [[11], |29|, 0, |30| are based on the chaotization 
of the intrinsic dynamics and they will not be repeated here. 
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7.2 Excitation and decay of a simple mode 

Here we compare the cross sections of various processes initiated in the channel c. They start 
with the excitation of the simple mode. The "elastic" scattering, c — > c, competes with the 
evaporation c — > e through numerous compound channels e. These branches are described by 
the amplitudes T cc and T ec , respectively, see eq. (|53|) . 

On the real energy axis, z = E, the uniform model leads to the inverse Green function 
Qoo(E) = F(E), compare eq.flD, 



'2(1 



x 2 y 2 ) 



At y = (no evaporation), the elastic cross section 



\T 



cc\2 



ll/\HE)\< 



7o 2 



E 



7o + T a y 



irE 



:i+x 2 



1 + x 2 y 2 



(62) 



ycot — ) ,,7/1 



-i 



(63) 



reveals fine structure fluctuations. In the case of small y ^ 0, these fluctuations are enhanced 
in a vicinity of the point E = eo due to the energy dependence of the imaginary part of 
F(E), eq.([62|). However, the fluctuations are washed away when evaporation becomes strong, 
■y ev ^> -D, so that y — > 1 and (|62"D simplifies to 



^(E) = E - e + - ( 7o + T s 



(64) 



Note that here the decay width 70 and the spreading width T s are combined into the total 
width of the resonance on the real energy axis. In Fig. 3 we illustrate the energy dependence 
of the elastic cross section a cc = \T CC \ 2 for different values of relevant parameters. 



Using the optical theorem, one obtains from eqs. 



-2ImT cc (E) = \T CC (E)[ 



and 



7e 



7o 



(65) 



for the total cross section initiated in the channel c. The fraction of \T CC \ 2 in the total cross 
section determines the branching ratio of the simple decay mode, 



B CC {E) 



7o 



lof(E) 



70 + j ev L(E) 7 /(£) + 7e 



f(E) 



(66) 



in agreement with the probabilistic interpretation of the function f(E). This function rather 
than its discrete counterpart fllBD is relevant when an actual reaction process is considered. 

The amplitude T ec (E), eq.(^), for evaporation in a given channel e after the simple state 
is excited in the entrance channel c, strongly fluctuates together with the exit amplitudes A e v . 
This amplitude vanishes in average. Assuming many uncorrelated statistically equivalent decay 
channels, we can use a natural statistical suggestion O, 



Taking into account eqs. 



(A^At) = S ee '5, ulev /k. 
we obtain 



1 2 7er 

&7, 



■L(E), 



(67) 



(68) 
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so that the corresponding branching ratio is equal to 



B' 



levL(E) 



''lev 


1 - f{E) 




lof(E) 4 


7et> 


1 - f(E) 



(69) 



Eqs. ( po|) and ( p9]) give £> cc +/c£> ec = 1 in accordance with the unitarity condition. The statistical 
ansatz (|67D is self-consistent because an equivalent approximation was in fact introduced earlier 
when the off-diagonal elements of the antihermitian operator w in the compound space were 
substituted by the average evaporation width, see the discussion after eq. ([TT]) . 

In the case of considerable evaporation and overlapping compound resonances, 'jev/D ^> 1, 
the branching ratios saturate at, see (|60|), 



B c 



7o 



7o + r s 



B ( 



k 7o + r s 



(70) 



For the saturated spreading width T s , these limiting values cease to be sensitive to the level 
density of compound states and depend on excitation energy or temperature only through the 
direct width 7 . Under such conditions, only the simple state with the total width 70 + T s , 
corresponding to the two possible ways of its decay, escape and internal dissipation, is seen in 
the scattering in the entrance channel c. Here again the background of compound states serves 
as a reservoir for irreversible decay, equivalent by its properties to decay into continuum. 



7.3 Reactions initiated in compound channels 

The processes started in the compound channels e, for example, driven by a nuclear interaction 
of heavy ions, can populate the simple mode through internal mixing. The corresponding 
amplitude T ce is the same as the amplitude T ec considered above, eqs. (|53"D and (^). The 
competing compound-compound processes are described by the set of the amplitudes 

A e ' A e A e '\/ V A e 

T e ' e (E) = E + E f^r^ E ^Ff- (71) 

v hi — e v v & - e v ^ & - 

The second term in ([fl]) accounts for the virtual excitation of the simple mode with the subse- 
quent deexcitation again via compound channels. 

To evaluate the total cross section of compound-compound reactions, 

a e _ £<| T e'f ), (72) 



we perform here the statistical averaging as in (p7[) . Neglecting the numerical corrections of 
the order 1/k, and using the notations of Appendix for the sums over the spectrum of the 
compound states, we obtain 

a e = -^{\S\ 2 + kS n + 2Re \g* 00 (V 2 )(SS 2 + kS l2 ) 



+ \Go \ 2 (V 2 ) 2 (\S 20 \ 2 + kS 2 n )}. (73) 



The terms proportional to the number k of open channels appear as a result of pairwise coherent 
averaging of random decay amplitudes. Only these terms survive in the limit fc > 1. Taking 
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the sums of Appendix in the overlapping limit y — ► 1 and recalling that in the same limit, 
according to (|68D, 

a ce = |T ce | 2 = ^oo| 2 7 r s , (74) 

we come, after many cancellations, to a simple expression for the total cross section of all 
reactions initiated in a generic compound channel e, 



<j + a 



rpce\2 _|_ V^ ^j ye'e|2^ 
1 



(75) 



2tt 

T 



(£ - eo) 



(70 + r s 



l27e?; 



This gives the branching ratio for the deexcitation into the channel c carrying the signature of 
the simple mode, 

B« = - ^ — . (76) 

The resonance at the simple state is suppressed by the inverse number j ev /D of the background 
states on the typical evaporation width. As it was discussed above, the observation of the signal 
of the simple mode in the reaction started in a compound channel becomes less probable with 
increasing j ev , in agreement with the kinetic arguments of pO|j . 

The same result can be expressed with the aid of the function /, eq. (|56|) . 



1 r 2 D 

^7or s |^oo| |T- 



(77) 



Here fa ce determines the fraction of the cross section of the process e — > c due to the intrinsic 
simple state; the denominator is the similar fraction of the compound-compound cross section 
due to the complicated intrinsic states, with no excitation of the simple mode. The right hand 
side of (|77|) is the resonance curve of the simple mode excited through the background (entrance 
factor r s ) and deexcited through its own exit channel (factor 7 ). The integral of the left hand 
side ratio over the energy region covered by the spreading width gives the inverse number of 
fine structure states in this region, D/Y s . 



7.4 Common decay channels 

One of the objections raised against the kinetic explanation [^] of the disappearance of the 
collective strength of the giant resonance is related to the possibility of preequilibrium excitation 
of the giant mode [35R. In this case the intrinsic evolution would start with the state which 



already carries some amount of collective strength. In our language such a possibility can be 
taken into account via the presence of the reaction channels a connected both to the simple 
mode and to the background states. For such channels, all amplitudes, Aq and A a v do not vanish; 
until now we assumed that, before the internal mixing, the simple state and the fine structure 
states have no common decay channels. For the case of the IAS, this situation is associated 



with the extrnal isospin mixing which is apparently of minor importance ||11|| . However, for the 
giant resonance this effect can change the situation. 

The common decay channels can be incorporated into theory without problems. Here we 
consider the simplest case of a single common channel which can be easily analyzed by the 
standard means. The corresponding real amplitudes will be denoted as ao and a v for the simple 
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state and bacground states, respectively. A many- channel case brings in many amplitudes of 
such type. Being uncorrelated, they should not lead to any effects of coherent enhancement. 

In the single channel case, the exact algebraic solution gives the matrix elements of the total 
Green function (B3T) 



£oo(£) 



E - e - £ 



V 2 



E 



— + 



a? 



QUE) 



E 



— + 



Guo(E) 
1 (i 



2 1 + (i/2)R° 



/oo 



E- 



~ 5 



E — 



2 1 



{i/2)R a 



(78) 
(79) 
(80) 



Here the renormalized amplitudes are introduced for the decay of the simple state through the 
channel a, 

«o(£) = ao + El^ (81) 

v & ~ e v 

and for the mixing between the simple state and the background including the intermediate 
continuum states, 



u v {E) = V v 



oc 



2 "l + (i/2)i2°' 



The analog of the i?-matrix, eq. (|35|) , for the a channel is 



R\E) = £ 



E-l 



52) 
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Using these exact expressions we evaluate the reaction amplitudes. We are interested in 
reactions starting in the channel a and ending either in the channel c specific for our signature 
of the simple mode or in any of the other channels, a or e. The elastic a — » a amplitude is 
given by (compare (]36f) 



R a 



+ 



'00- 



l + (i/2)R a \l + (i/2)R a / 
The deexcitation through the special channel c is governed by the amplitude 

G 



T ca = ^(Gooao + Gova v ) 



loGi 



00" 



^4) 



5) 



1 + (i/2)R a ' 

To make a conclusion of the importance of the excitation through the common channel, we 
assume that, similar to the amplitudes A e v1 eq. (|67|) , the new amplitudes a v are uncorrelated 
quantities with a large magnitude which contributes significantly to the total width 7 e „ of the 



compound states, (a ) 



la 



7 e „ ^> D. Using the estimates of the Appendix for the sums 



over the fine structure states in the overlapping limit, we obtain R a m —in^a/D). Thus, R a 
is a large imaginary quantity determined by the number of compound states in the interval 
7 ffl . According to the same estimates, the first term in ([84]) dominates, and the contribution of 
terms containing the sum with the cross products V u a u is relatively small. Finally, 



|T m | 2 ) 
|T aa | 2 ) 



7o7ol& 



00 1 



16) 
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The average partial width for the decay of the simple state into the channel a is equal to 



7o a EE(KI 2 > = ^ + ^r s . (87) 

''lev 

This, quite general, result shows that, in the case of the common channel capable of populating 
both, simple and compound states, the statistical branching for the deexcitation via the simple 
mode drops with the increasing level density p ~ l/D of compound states. Whence, strong 
common channels with 7 a ~ ^ ev cannot recover the disappearing simple mode. Weak channels, 
7 a ~ jev/k, are useless because of their small total cross sections. 



8 Conclusion 

In the paper we considered the most general properties of an open quantum system where a 
simple mode of excitation interacts with the background of very complicated states. Both, 
simple and compound, states are coupled to the continuum and have finite lifetimes. Internal 
dynamics (mixing) and external dynamics (decays) are intertangled in a nontrivial way. The 
intrinsic dynamics in the presence of the continuum are governed by the effective nonhermitian 
hamiltonian. The widths of intrinsic states modify the strength function of the simple mode. 
This view from "inside" has to be supplemented by that from "outside" for determining the 
observables measured in a real scattering experiment, such as cross sections and delay times. 
In a formal language, here we project the dynamics of quasistationary intrinsic states back 
to the real energy axis. The effective hamiltonian by its construction guarantees the correct 
properties of the scattering matrix including unitarity. Therefore it becomes possible to use 
the knowledge of internal dynamics in order to compare cross sections of competing processes. 

The general although schematic character of the analysis allows one to draw the conclusions 
concerning the manifestations of the simple mode in various situations. A typical example is 
given by the IAS which can be seen as a broad pole [|K]] or to be dissolved in the sea of the 



fine structure levels of another isospin. The analysis, analogous to that in |TT[ , confirms the old 



idea |ZT , |22|| of increasing isospin purity of the IAS at high excitation energy. The experimental 
data agree with this conclusion. The isospin purity is restored because of the very fast 
depopulation of the admixed background states when their decay width ^ ev increases compared 
to the spreading width T s of the simple state (IAS in this case). Here T s is assumed to be a 
slow changing or saturating function of excitation energy B as predicted by the analysis based 



on the chaotic character of the intrinsic dynamics p9 |. 

Such a consideration is not specific for the IAS and can be applied to other simple modes 
embedded into continuum. The giant dipole resonance is known to preserve its individuality 
up to high excitation energy or temperature []IJ. In particular, this is clearly seen in the 
observation of the nearly harmonic double-phonon excitations M. The new phenomenon of 
disappearance of the collective strength of the GDR [|IJ is still a debatable subject. Such a 
behavior was qualitatively explained in as a result of a shift of the kinetic equilibrium in 
favor of compound decays when the ratio jey/Tg increases. Our general quantum-mechanical 
analysis confirms this result. Moreover, we made the arguments which demonstrate that the 
conclusion is still valid when the simple mode can be excited from the reaction channels which 
are common for the simple mode and the background states. 

To complete the analysis, it would be interesting to consider the situation when several 
simple states can share the collective strength and the decay width into the channel which 
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signals the deexcitation of the simple mode. This is the case in the realistic calculation of 
the GDR. The collective peak accumulating a large part of the isovector dipole strength is 
shifted to high energies compared to the unperturbed shell model position. However, some 
strength is still concentrated at the unshifted energy. This "configuration splitting" leads to 
specific interference phenomena [16] which again can be described with the use of the effective 
nonhermitian hamiltonian. The distribution of the dipole strength and the width evolves with 
the increasing excitation energy which should be taken into account when the interplay of the 
internal interaction and external decays is considered. Typically, this results in the quenching 
of the collective strength and its redistribution in favor of the low energy component. These 
effects 



36 are discussed in 
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Appendix 

A uniform model of compound spectra 

Assuming the equidistant spectrum of unstable background states with the level spacing 
D and the decay width 7, which corresponds to -y ev of the main text, and substituting actual 
coupling matrix elements V 2 by their average (V 2 ), we have to deal with the sums as the trace 
of the Green function QTTD 

00 

S= E UE). (88) 



For the calculations of the scattering processes, the energy E is real and 

1 1 



The summation in (|88|) leads to 



E 



E-uD + (2/2)7' 



i9) 



7T 



S = — cot 



where the parameters are introduced 



x 



cot 



7T % 
— \E + -7 
D V 2 [ 



fuE 



7T X 



iy 



D 1 



ixy 



y 



tanh (Id) 



(90) 



(91) 



As the decay width 7 increases, the quantity y changes very rapidly from a small value 
y m 7T7/2D for isolated long-lived states, when 7/D <C 1, to a value exponentially close to 1 for 
overlapping levels, when 7/D 3> 1. In practice it is sufficient to consider just these two limiting 
cases. At small 7, the imaginary part of S is small, oc y ~ (Trj/2D), and the real part of S is 
equal to tcx/D as for stable levels |T7|]. In the opposite case of large 7/D, the real part vanishes 
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~ (1 — y 2 ) whereas ImS ~ —n/D. Both cases have a general meaning being not limited by 
restrictions of the uniform model. Thus, the result for the overlapping case follows immediately 



after substituting the ^2(E 



by the integral over the levels with a level density 1/D and 



using a small shift of energies into the complex plane. This expression is routinely used in 
statistical theory of nuclear reactions [11]. A similar consideration is valid for the sums as in 



taken at a fixed complex energy Ej = E — (i/2)Tj instead of running real energy 



( 13 ) and 
E. 

More complex sums can be easily analyzed in the same way. Here we give some examples 
used in the text (the notation S mn corresponds to m factors ( u and n factors Q so that the 
basic sum S = «Sio): 



S 
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V 



7T \ • 

d) 



(l-y 2 )(l + x 2 ) 
(1 + ixy) 2 

2tt y(l + x 2 ) 



(92) 



(93) 



<si 2 =£ic,rc 

V 

In the overlapping limit y - 



1 -y 2 ){l - x 2 y 2 + 2ixy) 



(94) 



v D'j 1 + x 2 y 2 

™ 1 +^ 2 _ _„.2\n ..2, .2 

D^i 1 + x 2 y 2 7 D 

1, these sums go to 0, 2n / ' and 2m /-D7 2 , respectively. The 
first sum ( p2|) vanishes in accordance with the fact that both poles in the equivalent integral 
are located on the same side of the real axis. The nonvanishing sums are proportional to the 
level density p = 1/D, i.e. they have a coherent component growing at high excitation energy. 
The sum S22 can be calculated as 

& = EICI 4 

V 

After simple algebra, we obtain 



d 2 



dEdE 



(95) 



E'=E 



4 T „ 2 „ dS 
S22 = — AmS + —Re— 



r 



dE 



(96) 



In the overlapping limit, the first term in ( p6[) gives Air/D^ 3 whereas the second one is propor- 
tional to (l—y 2 )/D 2r ) 2 and therefore it is small compared to the first term since the exponential 
smallness of (1 — y 2 ) over compensates an extra factor j/D. 
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Figure captions 

Figure 1. The strength function fi as a function of energy for the values of parameters 
T s /D = 100 and |7 — j ev \/D = 90. The solid curve gives the exact numerical solution of 
eq.(^), the dotted line corresponds to the approximation (p6|), the dash-dotted curves show 
the Breit-Wigner approximations (^9|) . 

Figure 2. The relative probability f(E), see eqs. (56-58), of excitation of a simple state 
through the channel c as a function of energy, E/D, and the evaporation width, ^ ev jD. The 
value of the spreading width, T s /D = 10 is chosen for illustrative purposes to make the oscil- 
lations along the energy axis clearly seen; the oscillations rapidly disappear as 7 e „ grows. 

Figure 3. Elastic cross section a cc in the channel c as a function of energy. The parameters 
in part (a) are 70/-D = 10 and T s /D = 5; the cross section is shown for different evaporation 
widths, •jev/D, which correspond to the values y = 0.3 (dots), 0.7 (solid curve) and 1 (dash- 
dotted curve). The situation with T s > 70 is shown, for the same values of y, in part (b) where 
Jo/D = 5 and T s /D = 10; note the different scale for the cross section. Part (c) shows the 
cross section for 70/-D = 100 and T s /D = 15 with y = 0.5 (oscillatory curve) and y = 1 (thick 
curve). 
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